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ABSTRACT 


\ 

This  report  contains  a Fortran  program  that  enables  the  user  to 
produce,  from  published  tables  of  normalized  element  values,  the  circuit 
diagram  and  amplitude  response  of  a desired  LC  ladder  filter.  Butterworth, 
Thomson  and  Chebyshev  characteristics  are  Included,  and  the  filter  may 
be  low-pass,  high-pass  or  bandpass  of  any  order  from  1 to  8.  Some  commonly 
used  sets  of  tables  are  Included  Internally.  Several  design  examples  are 
given. 
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INTRODUCTION 


Though  active  RC  networks  are  now  popular  in  filter  design,  there  Is 
still  a place  for  the  passive  LC  variety.  The  latter  is  Invaluable,  for 
example,  in  situations  requiring  a minimum  use  of  power,  such  as  an  installa- 
tion at  a remote  site  or  an  underwater  location  with  a long  cable  link  to 
shore . 


Tables  have  been  published  by  Weinberg  (Reference  1)  that  give  the 
normalized  element  values  for  Butterworth,  Chebyshev  and  Thomson  (Bessel) 
filters  of  various  orders.  This  information  is  sufficient  to  enable  the 
user  to  design  filters  with  little  difficulty.  However,  use  of  a computer 
can  further  reduce  the  amount  of  work  Involved;  to  a bare  minimum  in  fact. 

In  the  appendix  is  the  listing  of  a Fortran  IV  computer  program 
that  yields  circuit  diagram  and  plot  of  the  amplltude-vs-f requency  response 
of  a desired  filter.  The  required  input  data  are  entered  on  one  or  two 
DATA  cards. 

GENERAL  RF.MARKS 

I 

The  general  form  of  the  filters  considered  in  this  report  is  shown 
in  Figure  1.  This  network  may  be  a low-pass,  high-pass  or  bandpass  filter. 

It  can  have  any  number  of  reactive  elements  (order  of  filter)  up  to  eight. 

The  response  characteristic  is  selected  from  the  following  group: 

(a)  Butterworth  - maximally  flat  magnitude. 

(b)  Thomson  - maximally  flat  time  delay. 

(c)  Chebyshev  - equal-ripple  magnitude. 

The  Weinberg  tables  list  the  element  values  (in  units  of  ohms, 
henrys  and  farads)  for  low-pass  filters  that  are  normalized  with  respect  to 
frequency  and  impedance.  For  each  type  of  filter,  these  normalized  values 
are  grouped  in  sets,  each  under  the  heading  of  r,  which  represents  the  ratio 
of  input  and  load  resistances.  In  each  set  are  ten  rows  of  element  values, 
one  for  each  value  of  n,  the  order  of  the  filter.  The  maximum  value  of  n 
used  by  the  program  is  eight,  though  it  would  not  be  difficult  to  extend 
this  to  ten. 

The  basic  low-pass  filter  is  shown  in  Figure  2.  In  normalized 
form,  R^  » 1.0  ohm;  R may  be  0.0  ohm  (short  circuit)  or  a value  determined 
by  r.  If  n is  even,  r = odd,  r = filter  con- 

figuration may  be  as  shown  (n“8)  or  have  fewer  reactive  elements.  It  is 
important  to  note  that  the  labelling  order  of  the  element  values  taken  from 
the  Weinberg  tables  is  reversed  in  this  report.  In  a fifth-order  filter, 
for  example,  L5  becomes  LI.  This  should  cause  no  difficulty,  especially 
after  two  or  three  examples  have  been  studied. 

A few  words  of  explanation  are  required  with  regard  to  the  frequency 
normalization.  Both  Butterworth  and  Chebyshev  filters  are  normalized  to 
make  u>^  * 1.0  rps.  For  the  Thomson  filter,  however,  w * 1/t  , where  t 
is  the  time  delay  at  zero  frequency  (different  for  eacR  order°of  filter?. 
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As  will  be  seen,  the  computer  program  arranges  matters  so  that  all  of  the 
designed  filters  have  a specified  amplitude  response  at  the  design  cutoff 
frequency.  Usually  this  level  is  -3.0  dB. 

TRANSFER  FUNCTION 


The  transfer  function  of  the  general  network  shown  in  Figure  1 is 
easily  obtained: 

Let  E^  = 1.0,  R^  = 1.0 

h ■ '=2  ■ 

V = Z^I„+1.0 

c 7 2 
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g h 2 
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For  the  low-pass  filter  of  Figure  2,  the  Impedances  and  admittances  are  as 
follows : 
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Note  that  the  transfer  function  (1)  Is  applicable  to  filters  of  lesser 
orders  than  8.  For  a sixth-order  filter,  for  example,  Z_  " 0.0  and  Y * 0.0. 
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FREQUENCY  TRANSFORMATION  AND  FREQUENCY  SHIFT 

Transformation  of  the  frequency  variable  Is  used  to  convert  a low- 
pass  filter  into  a high-pass  or  bandpass  filter.  Obviously,  no  such  trans- 
formation is  necessary  if  a low-pass  design  is  required;  however,  in  some 
cases  an  extra  frequency  shift  is  essential,  as  will  be  explained.  Though 
the  actual  transformations  are  not  discussed  here  (see  Reference  1) 
the  consequent  steps  are  given. 

(1)  Low-Pass  Filter 


To  convert  from  the  normalized  filter  to  one  with  a cutoff  fre- 
quency (ij  and  a load  resistance  of  the  required  value,  proceed  as 
follows : 

(a)  Multiply  all  normalized  element  values  by  the  radian 
frequency  at  which  the  amplitude  response  is  down  3 dB. 

(This  step  is  necessary  because  of  the  way  in  which  the  Thomson 
filter  is  normalized.  For  the  Butterworth  filter,  = 1.0  rps 
For  the  Chebyshev  filter,  the  amplitude  at  the  corner  frequency 
is  down  an  amount  equal  to  the  declared  passband  ripple.  Thus, 
(jjj^  is  slightly  greater  than  1.0  rps  for  ripples  less  than  3 dB. 

A little  forethought  is  required  in  this  connection,  as  will  be 
explained  in  the  section  describing  the  program.) 

(b)  Divide  all  L's  and  C's  by  u . 

c 

(c)  Multiply  all  L's  and  divide  all  C's  by  K,  the  impedance-level 
parameter. 

(2)  High-pass  Filter. 

(a)  Obtain  the  normalized  low-pass  filter.  Including  multiplication 
by  03^. 

(b)  Replace  each  L with  a C equal  to  1/L  and  replace  each  C with 
an  L equal  to  1/C,  the  units  being  henrys  and  farads. 

(c)  Divide  every  new  L and  C by 

(d)  Multiply  all  new  L's  and  divide  all  new  C's  by  K,  the  impedance- 
level  parameter. 

(3)  Band-pass  Filter. 

(e)  Obtain  thr  normalized  low-pass  filter,  including  multiplication 
by  u,^. 

(b)  Divide  all  L's  and  C's  by  w ■ Uj-u-,  where  w.  and  are  the 
edge  frequencies  of  the  desfred  pass^band. 
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(c)  Multiply  all  L's  and  divide  all  C's  by  K,  the  Impedance-level 
parameter . 

2 r 

(d)  To  each  L add  in  series  a capacitance  C » 1/LWq,  where 

to  each  C add  in  parallel  an  Inductance  L = 1/Cw^. 

PROGRAM  DESCRIPTION 

It  is  not  intended  to  explain  the  program  line  by  line,  but 

rather  to  indicate  the  purpose  of  the  various  sections.  It  will  help  to 

know,  for  instance,  that  almost  half  of  the  main  program  and  all  of  the  sub- 
routines are  used  to  draw  the  circuit  diagram  and  label  it. 

Since  a knowledge  of  the  input  data  variables  is  important  to  an 
understanding  of  the  rest  of  the  program,  these  are  listed  immediately. 

PASS  - determines  the  kind  of  pass  band. 

1 - low-pass 

2 - high-pass 

3 - bandpass 

N - order  of  filter 

FC  - cutoff  frequency  of  low-pasa  or  high-pass  filter 

FI  - lower  corner  frequency 

F2  - higher  corner  frequency 

K - impedance-level  multiplier 

DF.BCOR  - decibel  response  at  cutoff  frequency 

SELECT  - chooses  element  values 

0 external  list  entered  by  means  of  DATA  cards 

1 Butterworth,  RA»0 

2 Butterworth,  RA*1.0 

3 Chebyshev,  RA»0  (1-dB  ripple) 

4 Thomson,  RA=0 

5 Thomson,  RA-1.0 

R - resistance  ratio,  r. 

Lj^,  C^,  Lj,  C^,  L^,  C^,  L^,  Cg  - element  values. 

These  variables,  corresponding  to  those  in  the  NAME  LIST  statement, 
are  entered  by  means  of  DATA  cards  placed  immediately  before  the  END  OF  DATA 
card.  The  list  is  largely  self-explanatory,  but  further  clarification 
will  be  given  during  the  remainder  of  the  program  description  and  the 
examples  that  follow. 

The  large  block  of  numbers  appearing  after  the  NAME  LIST  statement 
in  the  program  is  a three-dimensional  array  called  FILTER.  This  array 
contains  eight  rows  (n“l  to  8)  of  normalized  element  values  for  each  of 
five  sets  of  low-pass  filters,  making  forty  rows  in  all.  These  sets  have 
been  put  in  the  array  for  the  convenience  of  the  user.  Included  are  the 
five  cases  given  above  opposite  the  variable  SELECT.  If  SELECT  ■ 0, 
the  entire  block  is  skipped  and  the  external  normalized  element  values 
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selected;  If  SELECT  = 2 (for  example),  selection  is  made  from  FILTER  of  the 
row  of  normalized  element  values  that  will  produce  a Butterworth  filter  of 
order  N with  input  resistance  equal  to  load  resistance. 

There  should  be  no  difficulty  in  knowing  which  row  has  been  chosen. 

N is  easily  picked  out,  for  it  equals  the  number  of  non-zero  values  given 
in  the  row.  The  sets  are  listed  in  the  order  shown  opposite  SELECT.  As  a 
check,  the  selected  values  appear  in  the  print-out. 

These  normalized  element  values  are  used  as  a basis  for  the  design 
of  all  the  filters  dealt  with  in  this  report,  whether  they  are  low-pass, 
high-pass  or  bandpass  in  nature.  The  choice  of  this  option  is  indicated 
by  setting  PASS  = 1,  2 or  3,  respectively.  Other  required  information 
includes  FC,  the  cutoff  frequency  in  the  case  of  a low-pass  or  high-pass 
design,  or  FI  and  F2  if  the  design  is  to  be  bandpass.  The  impedance-level 
parameter  K is  also  necessary,  as  well  as  the  resistance  ratio  r,  which  is 
associated  with  the  chosen  element  values.  If  SELECTED,  r is  set  internally 
and  is  not  required;  nor  is  it  necessary  to  assign  element  values. 

The  remaining  piece  of  Information  is  DEBCOR,  the  decibel  response 

at  the  cutoff  frequency.  Although  this  is  usually  set  equal  to  -3  dB , there 

Is  no  reason  - with  one  precaution  - why  it  should  not  be  some  other  value. 

The  precaution  is  concerned  with  the  undulating  nature  of  the  amplitude 

response  in  the  pass  band  of  a Chebyshev  filter.  If  the  magnitude  of 

DEBCOR  does  not  exceed  the  ripple  magnitude,  an  erroneous  result  will  be 

obtained  for  w,  . 

k 

Given  a set  of  element  values,  the  computer  proceeds  by  a step-by- 
step  calculation  to  determine  the  amplitude  level  at  each  frequency.  When  the 
level  drops  below  DEBCOR  the  calculations  cease  and  the  radian  frequency 
at  that  point  is  designated  w^^.  This  result  can  be  negated,  if  desired, 
by  inserting  after  = w a card  reading  w^^  = 1.0.  There  may  be  occasions 
when  it  is  desirable  to  do  this,  but  for  most  designs  it  will  not  be 
necessary. 

Having  found  w.  and  knowing  w and  K,  the  program  then  calculates 
the  component  values  (R^s,  L's,  C's)  for  the  actual  filter  in  the  manner 
described  in  the  previous  section  of  the  report.  The  amplitude  response  of 
the  filter  is  now  calculated  and  the  results  plotted. 

The  considerable  length  of  that  section  of  the  program  dealing  with 
drawing  the  circuit  diagram  is  due  to  the  number  of  possibilities  that  exist. 
The  network  may  be  low-pass,  high-pass  or  bandpass;  it  may  be  of  any  order 
from  1 to  8,  odd  or  even;  the  input  resistance  may  be  present  or  not;  the 
element  values  may  be  given  in  henrys  or  millihenry s,  microfarads  or  pico- 
farads, with  any  of  a number  of  decimal  places. 

Despite  the  length  of  this  part  of  the  program,  inspection  will  show 
it  to  be  segmented  and  straightforward.  Once  the  pattern  is  recognized, 
understanding  comes  easily. 

Finally,  all  one  need  know  to  obtain  one  of  the  filter  designs 
; covered  in  this  report  is  what  Information  must  be  supplied  to  the  program 

J on  one  or  two  DATA  cards.  Some  examples  are  given  below. 
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Several  examples  are  given  to  clarify  the  use  of  the  program  and 
to  demonstrate  its  versatility. 

1.  Design  a sixth-order,  Butterworth  low-pass  filter  to  have  a cutoff  fre- 
quency of  9,500  Hz,  with  input  and  load  resistances  of  1,500  ohms  and 
6,000  ohms,  respectively. 

The  DATA  cards  should  read  as  follows: 

PASS=1,  N=6,  FC=9500.0,  K=6000.0,  DEBCOR=-3.0,  SELECT=0,  R=0.25,  L =0.0675, 
C2=3.1601,  L3=0.3130,  C4=6.4673,  L5=0.4567,  C6=6.3425,  *END*.  ^ 

The  element  values  are  obtained  from  the  Butterworth  tables  in 
Reference  1 under  the  heading  r = 1/4(=1, 500/6,000) . Note  that  the 
labelling  has  been  reversed.  Figure  3 gives  the  circuit  diagram  of  the 
filter  and  Figure  4 shows  the  amplitude-vs- frequency  response. 

2.  Design  a low-pass  filter  with  the  following  specifications: 

(a)  minimum  overshoot; 

(b)  cutoff  frequency  of  250  Hz  (-3  dB  point); 

(c)  amplitude  fall-off  to  reach  a maximum  rate  of  30  dB  per  octave. 

(d)  input  resistance  to  be  20,000  ohms; 

(e)  load  resistance  to  be  5,000  ohms. 

Enter  the  following  parameters  on  DATA  cards:  PASS=1,  N=5,  FC=250.0, 

K-5000.0,  DEBCOR=-3.0,  SELECT-0,  R=0.25,  Ll=0.2731,  02=0.0498,  L3=1.2690, 
C4-0.1084,  L5=2.8247,  *END* 

Note  that  the  specification  "minimum  overshoot"  dictates  the  use  of  the 
tables  for  maximally  flat  time  delay  (Thomson).  Since  n is  odd,  R = 
5,000/20,000  instead  of  the  reciprocal.  The  circuit  diagram  for  this 
design  is  given  in  Figure  5,  while  the  amplitude  response  is  shown  in 
Figure  6. 


3.  Design  a low-pass  filter  to  the  specifications  given  in  Example  2 except 
for  the  Interchange  of  input  and  load  resistances. 

Since  n is  odd  and  therefore  r * R^/R.  '^se  must  be  made  of  the 
reciprocity  theorem.  This  means  simply  that  the  design  proceeds  as  in 
Example  2,  assuming  the  input  and  load  resistances  are  as  given  there,  and 
then  the  network  is  turned  end  for  end.  Lift  the  bottom  of  the  5,000-ohm 
resistance,  making  this  the  input,  and  ground  the  (former)  input  side  of 
the  20,000-ohma  resistance  making  it  the  load  resistance. 

The  amplitude  response  should  be  the  same  as  in  Example  2,  except 
for  a much-reduced  flat  Insertion  loss. 
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4.  Design  a third-order,  Butterworth  higli-pass  filter  with  cutoff  frequency 
of  500  Hz,  input  and  load  resistances  to  be  12,000  ohms. 

Since  r=l,  the  Internal  element  values  can  be  used  by  letting 
SELECT=2.  The  DATA  card  reads  as  follows:  PASS=2,  N=3,  FC=500.0, 

K= 12000.0,  DEBC0R=-3.0,  SELECT=2,  *END* 

The  network  is  given  in  Figure  7;  the  amplitude  response  is  shown  in 
Figure  8. 

5.  Design  an  eighth-order,  Butterworth  bandpass  filter  with  corner  frequen- 
cies at  450  Hz  and  1,500  Hz.  To  minimize  insertion  loss,  the  input  resis- 
tance should  be  omitted  (RA=0,  assuming  a source  impedance  of  zero  ohms). 
The  load  resistance  is  8,200  ohms. 

Input  the  following  data:  PASS=3,  N=8,  Fl=450.0,  F2=1500.0,  K=8200.0 

DEBCOR=-3.0,  SELECT=1,  *END* 

The  network  is  shown  in  Figure  9;  the  amplitude  response  in  Figure  10. 

6.  Change  the  previous  design  to  that  of  a Chebyshev  filter  with  1-dB  ripple. 
All  other  specifications  remain  the  same. 

Simply  change  to  SELECT=3.  Figure  11  shows  the  circuit  and  Figure  12  the 
amplitude  response. 


CONCLUSION 

A Fortran  program  has  been  provided  as  an  aid  to  the  design  of 
LC  filters.  Its  use  in  conjunction  with  tables  of  element  values  given  in 
Reference  1 facilitates  quick  designs  of  Butterworth,  Thomson  and  Chebyshev 
filters  of  various  orders  and  pass  bands.  For  each  design,  a circuit  dia- 
gram is  drawn  and  the  corre.spondlng  amplitude  response  curve  is  plotted.  The 
entire  process  can  be  completed  in  minutes. 
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Figure  8.  Amplitude  response  of  filter  shown  in  Figure 
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Figure  10.  Amplitude  response  of  filter  shown  in  Figure 
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Figure  12.  Amplitude  response  of  filter  shown  in  Figure  11. 
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0*0000. 

0 

* 0000/ 

0 

ocoo. 

0 

0000. 

>0 

.6180/ 

1 .61 8o / 

?.0000. 

1 

.61 80. 

0*6180/ 

0 

*0000/ 

0 

OOOO. 

0 

OOOO. 

>0 

.5176/ 

1.4142/ 

1.9319/ 

1 

.9319. 

1*4142/ 

0 

*51 76 / 

0 

OOOO. 

0 

OOOO. 

>0 

.4450/ 

1 .?470/ 

1 .8019/ 

? 

.onoo. 

1*8019/ 

1 

*?47n/ 

0 

4450. 

0 

OOOO. 

>0 

.390?/ 

1»11  11. 

1.66?9/ 

1 

.9616. 

1*9616/ 

1 

*6629/ 

1 

1111. 

0 

3202. 

>0 

.5088/ 

0.0000. 

0.0000. 

0 

.0000. 

0*0000/ 

0 

*0000/ 

0 

ocoo. 

0 

OOOO. 

>0 

.9957/ 

0.91 lO. 

0.0000. 

0 

•0000. 

0*0000/ 

0 

*0000/ 

0 

ocoo. 

0 

OOOO. 

>1 

.5088/ 

1.333?. 

1.0118. 

0 

•OOOO. 

0*0000. 

0 

*0000/ 

0 

ocoo. 

0 

OOOO. 

>1 

.3817/ 

1 .9093. 

1 .41?6. 

1 

.0495. 

0*0000. 

0 

•0000/ 

0 

OOOO. 

0 

0000/ 

>1 

.6652/ 

1 . 59q8 . 

1.9938/ 

1 

• 4441. 

1 *0674. 

n 

*0000. 

0 

OOOO. 

0 

OOOO. 

>1 

.3457/ 

?.049l / 

1.6507/ 

? 

•0?70. 

1 * 4601 . 

1 

*0773. 

0 

0000. 

0 

ocoo. 

>1 

.7118/ 

1 .64  85. 

?.119?/ 

1 

.6736. 

?*0**3T/ 

1 

*4694. 

1 

o83?. 

0 

OOOO. 

>1 

.3691/ 

?.09?2/ 

1.7021. 

? 

•1453. 

1 *6850. 

2 

*0537/ 

1 

4751. 

1 

0872. 

>1 

.0000/ 

o.cnoo. 

0.0000. 

0 

•OOCO. 

0*0000. 

0 

*0000. 

0 

OCOO. 

0 

OOOO. 

>1 

.0000. 

0.3333. 

0.0000. 

0 

•OOOO. 

0*0000. 

0 

*0000. 

0 

OOOO. 

0 

OOOO. 

>0 

.8333/ 

0.4800. 

1.6667. 

0 

•OOCO. 

0*0000. 

0 

*0000. 

0 

OOOO. 

0 

OOOO. 

>0 

.7101. 

0.46?7. 

0.2899. 

0 

•1000. 

0*0000. 

0 

*0000. 

0 

OOOO. 

0 

OOOO. 

>0 

.6?3l  / 

O.'tPlS. 

0*3103. 

0 

. 1 948. 

0*0667/ 

0 

*0000. 

0 

OOOO. 

0 

OOOO. 

>0 

.5595. 

0.3821 / 

0.3005/ 

0 

.2246. 

0*1400. 

0 

*0476/ 

0 

ocoo. 

0 

OOOO. 

JTSCEDDC  PAGE  BLANK.NOT  FILMED 
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>0  *51 1 1 i 

0.3487# 

0.ER27# 

0.2288# 

0* 1 7o4# 

0 * 1 05S# 

0*  03S7# 

0*0000# 

>n*^T32^ 

0.3212# 

0*2639# 

0*2227# 

0* 1 8n6# 

0*1338# 

0*o823# 

0*0278# 

>?.0000j 

O.QOOO# 

C.OOOOi 

o.nooo# 

0»0000# 

0*0000# 

0*0000# 

0*0000# 

>0 • 4E26* 

1 .5^74# 

0.0000# 

o.noco# 

0*0000# 

0*0000# 

0*0000# 

0*0000# 

>n.l922i 

0 . 552R  # 

1 .2550i 

0*nn00i 

0*0000* 

0*0000# 

0*0000# 

0*0000# 

>0»1104j 

0.3181# 

0.5116# 

1 .n598# 

0*0000# 

0*0000# 

0*0000# 

0*0000# 

>n»0T18j 

o.2n9o* 

0.3312# 

0 * 4R77  # 

0*9303# 

0*0000# 

0*0000# 

0*0000# 

>0»0505« 

0. 148n# 

0.2364# 

0*3158# 

0*^116# 

0*8377# 

0*0000# 

0*0000# 

>0»0375« 

0*1104# 

0*1778# 

0*2378# 

0*2944# 

0 * 3774# 

0*7677# 

0*0000# 

>n • 0?89# 

0.0*5S# 

0*1387# 

0 * 1 267# 

0*2297# 

0*273S# 

0*3446# 

0*7l 2S/ 

L1 ■FILTER( l^N, SELECT) 

C?«ElLTEPf?,N, SELECT  ) 

L3-ElLTER(3iN, SELECT) 

C4»eilter(l#n, select  ) 

LS»EIlTER(S#N, SELECT) 

CA-EIlTER(6,N. select  ) 

L7-ElLTERf 7, N, select  ) 

C«-EIl.TER(8/N#SELECT  ) 

inn  c^intimjE 

wRIEEm08i91  )PASS,N»EC/Ei,r?iKiDEBCPP/SEl  EcT/P 
91  EBR'IATfSiA/SEl  n*2i  Ift^E7*3) 

write  (1  C8>9?  )L1  #C?#L3/C4,i.S,C6/L7#C« 

9?  E^R“1AT  <8E8.4  ) 

EInD  w< 

Or*  1 • 0 

P«N/2 

E«2*P 

rE(E.Ef:•^.9'?.P.EQ.O•0)C,P  Tft  16 
ra-pr/r 

r,9  T9  ?9 

1h  ra»rb*r 

29  CRNfiMjE 

EREQLC,  ( 1 ) .-P.7981  S 

D9  T.li160l 

erE'iji  I )»in.n**EREoLri(  I ) 

W«6*?832*EREa( T ) 

S.CMPlv (O.C/W) 

71 -LI *S 

v?»C2*S 

73»L3*S 

V4»C4*S 

7S»L5*S 

Vft»C6*S 

77»L7*S 

V8«C8*S 

I?.yg+1 .o/RP 

vc»Z7* i?+i #n 

IH-Y6*VC 

lr,«lH+r2 

VR-Z5* IG4VC 

Ir»v4**vB 

rr» ie4 iG 

VA«73»  TE-t-VR 


X 
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ro»  vA 
1 1 ■ ID*IE 

FI ■ (Ra+ZI )*ri+VA 

G(  I )»CAHS( (OA  + RR)/(RB*Fl ) » 

DFBEE  ( T ) ■?0.0*AL'’G10(Gf  I > ) 
TF(CERFE(  n .LE.DEPC9R)  r,0  TP  5 
FREQLGf  Ul  )«F»FOLn(  I )+n«OnPB 
h CPNTInljE 
«S  CPNTImjE 

WK»^ 

PBTArg  actual  C‘'Mppnfnt  valufs 
El  ( 1 ) -I  1 J EL  (?  ) «c? 

FL(3)»l3j  EL(^i»C4 
F|  5;  ELf6)«CA 

Fl  (7). I 7,  EL (8<-CP 
Dp  6 

Fl  F*^(  J)-O.C 
b CPNTiNtjE 
wr»P*?83P*FC 
V1«AL0G1O(FC) 

V?«AInT( AL9Rln<FCl  1 
Vl -vP 

lr( v3.l E.0.S)FREQI  G<1 )«V?*?.C 
fR(  V3.r,T.n.B)FREQLQ<l)»V?-1  .0 
Ra»RA*K»  rr»rr#k 
If(  .NBT.  (PASS.FQ.pn  GS  TP  R 
OP  7 

rF(EL(?*J-n  tFD.0.0)G9  Te  R 
Fl  )■!  .0/rEL(2*J-l  )»wit*wc*l<) 

rr(EL(?*j»  tFO.OtOiQ9  T0  7 
Fl  (?*j)»K/(FL( ?*J^*WK#wC) 

7 CPNTInuE 
GP  TO  1 P 

9 TF(RASS»EQ. 1 ) GP  T0  10 
wr"6»?R3P#(F2-Ft ) 
wnS0»3R»P78P#F1#F? 

Fn«SQRT (Fl *F?) 

V1 •ALenlOfFO) 

V?»AINT( AL9GlO(FCi » 

V3»Vl-vP 

ir(V3,lE.o«5)FRFni.G(l  )»V?-?,C 
rF(V3.GT,o.B)FREOLG(l .0 
in  DP  11  .i^liP 

FL (?*J-1 )-EL (?*J-1 )*K#W</wr 
Fl  (2»J^«EL(?»J>»Wk'/(*<*WCJ 
11  CPNTImjE 

TF(PASS»F0.1 ) G9  T9  IP 
OP  1? 

IF(EL(.|)  .Fa.O.niGP  T9  13 
Fl  ( J*«>"1 .0/(FI  (J)#WOSO) 

GP  T9  1 2 
13  Fl  ( J+8  > ■0.0 
1?  (!PNTINUE 
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PalCUlA'F  and  pi  8T  amplitude:  PPSP8NSE: 

1 A 15  T » 1 . API 

F«ED(  n«ln*0**PPEpLG(  1) 

A(.6.?a3?*PRFQ(  T ) 

S.CMPt  y(C.C.w) 
r,fl  T8  (1/?>3)PASS 
1 DH  51  ,j«l,a 
7PL  ( J ) .F-L  ( J 1 *S 
SI  rsN'^lMlE 
D'l  ■'■P  l9 
? Oft  17  j«i 

Ir(FL(,l)  .FQ.O.niGP  T0  17 
7Fl(j)  »1  •n/(E[.  (J)  *S) 

17  rsNTlMijE 

G0  79  ^9 

? OB  18  j>1/8 

7FL ( J) «Fl ( J1 * ( P+wpSO/S) 

18  DbNTImuF 

19  !?»ZFl,  ( 8 ) + 1 .0/Pn 
Vr-ZEl  (7),lP-n  ,D 
Th»ZEl  f 6 ) *VD 

Tr,-  IM*I8 

vr«zEl (5 1 * iR+vr 
1F»  ZEl  f ^ 1 *VR 
rr-IF^iG 
va»7El  ^ 3 ) * if ♦v/r 
Td»ZEl ( ? ) *va 
1 1 • ID^IF 

F1 • ( 9a^ZEl  M ) ) * I 1 +VA 
G( I ) •CABSf 1 .C/Fl  ) 

DFBEF( J ) •?C.0#ALPplC(3f I n 

I 7E  ( 1 C8,  Ep  IFPFOi  G ( I ) iFRFO  ( I ) / 0(  I ) iO'^BFF  (1  1 
?n  FbR>^AT  (aF19«A) 

IF(DEPFE( I ) .LF.-6p.C)DF3EFf I).-60»n 

FPEGLr,  ( Ui  ) .fpftlg(  I )>n*ri 

IS  CPNTImjE 

FrEOLG ( AC? ) -FRFDLG ( 1 ) i FRFDLG( 403 i «p. a 
FRE0(An2).ir.tc**F«FDLG(4n?ii  frfoiaogi.p.a 
DFBEF  ( aO?  ) --6P.C  j i^EbEE  ( APR  ) .1  0#0 
Tall  LlNFfFREDI  G/DEBEE #401, 1/0/0) 

TALL  AylSrO./O./' amplitude  (DB)  '/14/7.,pp*,DErFF(AP?)/DERFE(4p8) i 
CaLl.  LGAXS(P*/P./  'FREQUENCy  (H7  ) ' / -1  4/ 1 n.O/O.P/FREo  ( 40?  ) / 

>FRE0 ( 4o3  ) ) 

OR  ?1  .)■  1 , 4 

irt.NBT. (PASS. rO.?))ELFMfp*, 1-1). FL<?#j-1i 
IF-(.N97.(PASS.F0.?))ELEM(?*j).EL(!?*J)*1Pi.0**A.P 
lF(PASS.En,?)FlEM(?*J-1)«Fl  (2*J-l)*10.P**6.p 
IF(PASS.Eo.P)F|  EM(?*J).EL(?*J) 

IF(PASS.F0,3)EI  EM(2*jf7)»Fl  ( ?* J + 7 ) * 1 0 .n**6 .p 
IF(PASS.EOO)EI  EMf?*JAfi).Fl  (2*Ja8> 

?1  OPNTIm)E 

write  M 08/??)RA 
??  FRRMAT(F1o,1 ) 

OR  ?3  .i-l.a 

WRITE(108/?4)FI  EM(j)/ELEM(,u8) 


k 
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?4  F«R‘^AT(2Fi?.6) 

?3  r^NTlMlE. 

-<iWITEmO«,2'^)RR 
pe^  FhR'^AT  (Fl  C.  1 1 

r‘^R“AT(5FiO.?) 

DRA^  CIRrulT  DIAGRAf- 

r.ALL  PI  f3T  (1  4,n»  A.0<-3l 

x.0*0;  v.ntCj  CAt  I 

call  CTRCl  ( X» Y/C#0* 360«0intC5i0»CSi0»m 

TF(W.3T.0)G')  7P  ?7 

x«x*l.ni  CALL  PLBT(X#Y«?) 

Ir(PASS*Fn.3)A.-l  .0 

(PASS.ro.?)  )A«-c.A 
Call  P1.3T  , a»  o.(^* -p  ) 

GH  TB  P8 

P7  CBNTIniuE 

x»x>o.a;  cam  plbT(X»y«P) 

Call  SRFSH(yjY) 

x.x+o.P)  Cam  plbt(X,y#?) 

RA<«Ra/100C.O 

V«AINT(AlPG10(RA)  ijU-'s-v 

Call  m;mrfr  { .7*;^n.o7>n*lr«RA</o.CiU) 

CALL  sy"^rbl  ( 1 .?> n.c7/o • In# ' < ' / o*Oj  1 ) 

IF ( PASS. Fq. ? )R. -0.4 
iF(.NBT.(PASS.Fn.?))B«n.n 
CALL  PI  'll  ( Ri  O.n#  -?  ) 

?«  CBNTlMiE 

r,B  TO  ( icno#  Pnno# 3C00 )P*ss 
icon  DPI  1001  j-i/P 

i-aLl  pi  BT  m .6#  n.Oi  -3  ) 

X.O'Oj  X.fi.nj  call  PLBT(X#Yi3) 

M.AlM(ALPG10(10000.*FLF^'rP*J-n)) 
T.AINT(AlPG10{ 10Cn0000,*ELFM(2*j) ) ) 
fall  si  C 
icri  cpM'ImjF 

IF («.EO«n.O)LFN. ?.?♦!. fe*P 
Ir{P.GT.0.0)LF\»p.84l .6*P 
If(E.EO*N)G«  TP  1 POO 
CALL  PI  HT ( 1 .6#  n.C# -3) 
x.0*0i  >'»n.n;  CALI  plpt(x#v#3) 

‘^.AlNT(ALBG10(10CnO.*ElEM(?*P+l))) 

call  si 

IF(R.FO.O.O)LFN;«?.?4.1.6*(P+1  ) 

IF  < R.GT .0.0 )LFN»?.«  + 1 .ft* (P+1  ) 
iPnn  CPM^lMjE 

GP  ■'■p  400n 
Ponn  DP  ?ooi  j.l#P 

CALL  PL9T ( 1 .6#n.O* -3) 

X.O.Oj  Y.n.Oi  CALL  PLRT(X,Y,3) 

L ■ A INT  ( AlmGI  Of  1 OCnOOOO.#El.FM(P*j-l  ) ) ) 
0«AINT(AlP610(1 OOoO*  *ElEM( ?* j ) ) ) 

CALL  SCL 
P001  CPNI’INuE 

lF(P.En.0.0)LFN«?.?+l.A*P 
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ir  ( =!.GT.0.0  )LFN«?.a*l  «6*p 
Ir(E»FO#N)Gft  TH  ?P00 
CALL  PLSTM  .fe#n»C.-3) 

x.O'Oj  v.n.nj  call  plbto<»V/3) 

w. aint(Albgiomocooooo.*fi  rM(?#P4.l)  n 
CALL  sc 

If(P»FO*C.O)LF\»2.?4.1«A*<D4.1  ) 
lF(P»GTtO.CiLFN«?.S*1.6*(P+1  ) 

??nn  CbMInuE 

GB  TB  400n 
soon  DB  3001  J«l»p 

CALL  PI.9Tf2.0*n.0»-3) 

x. 0»0j  Y.OtOj  call  PLBTfXiY/3) 

M.A  INT  ( AlBGI  0(  1 00f'0**ELEM(  ?*j-1  ) ) ) 

L-AINTj ALBG101 10Cn0000.*FLFM(2#j+7j ) ) 
Q.AINTfALBGlOf 1 0000 • *ElEm f ?• J + 8 ) ) > 
T,AINTfA(_fl(5iO(  1 0000000.  #El  FM  ( 2*  J)  ) 1 
Call  SLCR1 

3001  CbNTInuE 

fF ( P.E0»0.CiLF\.?.6+2.0*P 

TF ( R.RT.0.0>LFN»3.2+2»n*P 

TF(E.En«N)GA  TB  3?00 

CALL  PLOT(2.0»0*C*-3) 

y»0«0i  Y.o.nj  CALL  PLBT(XiVi3) 

M.A  INT  f ALBGI  OM  00o0.#EL  EM(  ?*P+1  ) ) ) 

L-AIMT  ( AlBGI  0(  1 0000000.*ElFM(2*P4.9  1 ) 1 

Call  slcbp 

lr(P.En.0.0iLEN»P.2*2«0*(P4.1  ) 
rF(P.GT.0.0lLFN»2.g42.0*(P4.1  ) 

3200  CBNTIMjE 
4000  CBNTImjE 

IF(PASS.EQ.3»AND.F.EQ.N)GB  T6  30 

x.x+0.6)  Cali  plbT(X,Yj3> 

V.Y»0.65j  CAlj  PL0T(X,Y#2) 

CALL  SRESv(X#Y1 

V.Y-0.65i  cam  PLBT(X,Y,2< 

X.X-O.AJ  Y-V  + 2.U  cam  PL©T{X',Y«3< 

X.X*!.?;  CALI  PLBTlXfYsP) 

GB  TB  35 

30  continue 

X. X+l.Oi  CALI  PLBT(X,Y«3> 

Y. Y-0.ft5j  CAM  PLBTfX«Y,2l 

call  SPESv(X/V1 

Y.Y-0.A5)  CALI  PLBT(X#Yj2i 

X.X-1.01  Y»v^?,ij  CALI  PLBT(XfYi3l 

X.X*!.#.;  CAlj  PL0T(X,Vi3^ 

3S  continue 

Call  CI«CL  fX.V/l«0*0A-18o.O#O.05#O»o5»o.0> 
V.Y-2.1J  cam  PLBT(X,Y,3i 

CALL  CIf?CL  J X#  V#  lA0t0#-18O.0#0.05#n.O5,o.0l 
X.X-LEN  J call  PLBT(XiY«2l 

CALL  CIRCL  ^X/Y/0.0a360»0*  O»05/Ot0*5<  0.0) 
RRK-RB/1000.0 
V-AINT(ALOGIOIRR) ) J U»5»V 
RRXl  liX^LEN-OtA? 


Pr)(2«RRX1 

Call  MJf^RFP  tRRXl , -1  • 0«  1 *RRi<#0.0#U» 

Call  Sv'^RftL  (RRX?, -1  • 1#  0. 1 * 'K  • iO.O/ 1 ) 

re  { 3i , 3?i  )pass 

31  0.1.0*AlM(AL*lrilCfrC)  ) 

CALL  P(  »Tf x#-?,«^-3) 

x.c*0i  V.n.ni  CALL  PLi^TfX/V/P) 

x»  < LE\/2 •r ) - 1 . nR 

Call  5v'^RRl  (x«C»Pa  •LRw-PASS  F T L TER  • , 0 . 0/ 1 R ^ 
X.  ( lEN/?  . 0 ) - 1 • n*^  J v»»0 ‘O'?  J call  PL'iT  ( y < V#  T ) 
X«X>?.1J  TALL  PLBTfyiV,?) 

Cal^  sv^^r^l  fO.p^c.R^O'i/  'rr.'*n«nA3i 
Call  MjmrFR  (r.R#n.5/0.lAPr»c»0A-l  1 
Call  SY“ReL  ( .R5^n.i*DjO»Rjn*ii 'RZ'/n»n»p) 

GO  ■’’9  3“ 

3?  n.l.OvAlNTI AL9g1C<FC)  ) 

CALL  PL^T  (Xf?.R/-3) 

X.O'Oj  V.ntOi  CALL  PL9T(X#Y#3) 

X«(LEN/2»n)«l»1? 

CALL  Sy^HSL  (X*C»Oa «1A# 'HIrm-PASS  F I L TFR ' a n . Oj 1 A 1 
x«(LE\/2.0)-l.l2)Y»»C«0'ij  CALL  OL.9TfyAV,G) 
X»X>?.?A]  CALL  PL9T(y»y,p) 

Call  Symp«»l(0«PaC.5a0»1a  'Pr»'A0*0A3) 
call  MiMRFR  (n.Rin.®!AO.lArr/C.O#-l  1 

CALL  Sy“R9L  ( n»li  *MZ*  #n»n»P  ) 

GO  T9  GL 

33  n.l  .0-*-AlM(  AL9G10(F1  ) ) 

CAI.L  PI  9T  f Xa  -P.Pa  -3  ) 


34 


X«0*0j  Y.n.Oi  CALL  PL9T(XaVa3> 

n-l*0>AlNT(AL9Gir!fF?)  ) 

X.  (wE\i/2.0)»l  .HR 

CALL  SYMR9L  ( Xa  0»Oa  • 1 4>  ' RAK/DPASR  F T lTER  • a C . Oa  1 S ) 
X.  ( LE\/2#0  > -1 . n*?)  V»-0»05  j CALL  PL9T(xaYa3) 
x«X*2#i;  CALL  PLOTrxAYAP) 

Call  symrol  (nA?An»5A0.lA »Fi ■»a0.0a3i 

(0ASA0»5A0»tAF1 aO»Oa»1 1 
( AR54-n#l*DA0ARA0»lA  'HZ*  AO»nAP) 
(0.2An«35A0AlA*F2«'A0A0A3) 

(G.Ra .35a0.1aF?aC»0a»1 1 
( .R5-4-n»  1*Da  0»3Ra0»1  A ' H7  * A C.Oa  2 < 


MJ^IBFR 
SY*^B9L 
RYMB9L 
number 
SyMRBL 
continue 

CALL  NwPLOT 
STOP 

fnd 


CALL 

CALL 

•'ALL 

CALL 

CALL 
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9nl 


1 


9U0«eUTlNF  «;RE:RH(y,Y) 
901  M0P-1  <0 
y«x+0.n*» 

Y.Y-fO  .n^ 

Call  Pi.eT(x»Y#p) 

x»x+o*n’i 

Y«Y-o.n« 

CALL  PL9T(XjY,?) 
x»x-t-o.n^ 

Y.Y  + O .o'* 

CALL  PL0T(XiY#?) 

Continue 

pfTuRn 

fmD 


SuBRftUTiNF  SRESV<XjY) 
OR  1 I.liB 
X«X*0tO^J  Y«Y-n*CAl 
X-X-0*n8j  Y«Y-n»Onj 
x«x*o.n'»j  Y.Y-n.C4j 
continue 

PFTURN 

FnD 


CAI  I PLRT(  X,  y,?\ 
CAl.i  Pl6T(X*Y*?» 
CAI  I PLRTfXf  Y#?l 


SUBR0UTl\F  SLHfX/V) 

CALL  CTRCL  f X« Y# lfl0.0^-4B.n» 0» 1 »0» 1 /0«n 1 
X.X4.0«1707i  Y.Y-0.0707iCA|,|  PLRTjX,Yj3» 
CALL  CIRCL  f Xi Yi??5.0/-^5.n/0»l/0»i /0»0» 
X«X-*-0.1  ^14  J CAu  PL0T(XjY«  3) 

CALL  circl  <XiYj?p5.0/-45.n/C«i>o»i/n.n) 

X«X*0.1414j  CAl.1  PL0TfXjY#3» 

CALL  CIRCL  f X# Y/ ??5«0/-4B.Oi C» 1 »0» 1 /n»n ) 
X-X+0.141AJ  CALI  PLST|X»Y#3> 

CALL  CIRCL  (XiY#?P5.0/0«0«0. 1/0. 1/0.0) 
X.X+0.1 707i  Y»Y*0.0707|CAlI  PL9T(X/Y/3) 

Return 
end 


SUB^PUTINF  FLVfX/Y) 

CALL  CIRCL  f X/ Y/9o,0/ -1 3R,o^ 0.1 /O. 1 /O.O) 

X. X-0.0707J  Y.v-0.1707jCA|  I PLBTIX/Y/3) 
CALL  CIRCL  f X/Y/ l35.0/-13B.O/0.1/0.1/0.n> 

Y. Y-0.1414J  CALI  Pl«T(X/Y/3) 

CALL  CIRCL  ^x/Y/l35.0/-13B.^/0.1/0.1/0.n^ 
Y.Y-0.1414)  cam  PLeTlX/Y/3) 

CALL  CIRCL  fX/V/135.0/-13B.n/0.1/0.1/0.n< 
Y.Y-0. 14141  cam  PL0T(X/Y/3) 

CALL  CIRCL  fX/Y/135.0/-90.O/C. 1/0*1 /O.O) 
X»X*0.0707I  Y.Y-0.17071  CAM  PL0TIX/Y/3) 
RFTURN 

FND 


- A9 


SURRSUTINP  SL 

X,YiM,L*QiT#J#EL  EMf  if,) 
y»X-*-0.'‘H7?  j CAlI  PL9TfX,Y>?» 

TALL  SI H(Y/Y) 

y»X+0.4l7?i  CAl  I PL0T(XiY«2j 

RR  T9  f l<?>3<4#f5)M 

1 CALL  NUMRERf  0.R5#0»H*O«1.n*  1000»*ELEM(?*j-l  i<n.n*  n 
CALL  SYMBfiL  (n#90j0.14i0.in*  »MH**0t0>2>  I Qft  Tfl  f, 

2 CALL  MUMRER(0.B0#Otl4#o«lO#1000t#ELEM(?#,rt  wO.Oill 
CALL  SYMBOL  (n.95jO»U#0»1n* 'MH'iO.O#?)  I GB  TB  6 

3 tall  NlJMRERf  0«B5>0»1^aO»10*  lOOO»#ELEM(?*.rl  ) 

Call  symbol  {n»9o/o»l^* o»in* 'mh«#o»o*2) i go  to  6 

4 CALL  NUMBER(0«R5<O»l^AO»in»FLEM(2*.j-l  )#n.Oj2) 

CALL  Symbol  <i.00iO»l4*0»in/ ) i rO  to  a 

5 CALL  MjMRER  » 0 • B5#  0 • 1 0 • 1 0#  FLEM  ( 2#.)«1  ) # n • 0<  1 ) 

Call  symbol  n.oo*o»i4,c»ini 'hi#o«o<i  ) 
b continue 

x«l*6i  Y.n.nj  CALL  PL0T(X/Y»3) 

rfturn 

END 


SuRROUTlNr  SLC 

common  X,YiMiLiG<T/J#ELEMM6) 

Integer  t 
Call  sl 


Y-Y-1 

.0) 

CALI 

PLOTiX.Y.2) 

X«X-0 

• 1 5 j 

cali 

PLOTi X, Y/ 3» 

X«X>0 

.31 

CALI 

PL0T1X,Y/21 

Y.Y-0 

• 1 1 

cali 

PL0TiX/Y/3i 

X.X-0 

• 3 1 

CALI 

PLOTiX.Y/2) 

x»x+o 

» i 5 j 

CALI 

PL0TiX,Y/3> 

Y.Y-1 

.n; 

CAu 

PLOTl X.Y/21 

Y«  Y-»-2 

. 1 1 

CAl  1 

PL0TiX/Y.3> 

IF(T. 

LF«4)ELEMf 2*J) 

• ELFM(?#,i)  »10.0##6.0 

GO  TO 

11/2/3*4. 5a6>7#8)T 

1 

CALL 

number 

11.25/ 

-1 

.025/0.1 

/ ElE*^  1 2#  J 1 / 0. » • 1 ) 

CALL 

Symbol 

11.20/ 

-1 

• 1 75a  0. 1 

/ *PF'a0.a2)  iRO  to  9 

2 

CALL 

number 

1 1 .20A 

-1 

.025a  0. 1 

/ ElEm  ( 2#  J ) / («|,  , ^ 

CALL 

SYMBOL 

1 1 . ?0/ 

-1 

.175A0.1 

/ *PF'/0./2)  iRA  TO  9 

3 

call 

number 

11.10/ 

-1 

.025  a 0.1 

/ELEM(2«j)/n./*l 1 

CALL 

SYMBOL 

1 1 . 1 5/ 

-1 

. 1 75  a 0 . 1 

/ 'PF«/0.a2) |G0  TO  9 

4 

CALL 

number 

( 1 .00* 

-1 

•025*0 . 1 

/ElEH(  2*J>  /fi, . »i  ) 

CALL 

symbol 

11.10/ 

-1 

.175A0.1 

/ »PF • / 0./2 ) 1 rO  to  q 

B 

CALL 

Number 

1 O.ROa 

-1 

. 025 A 0.1 

/ ElEm( 2* J ) / R. , 4 ) 

CALL 

symbol 

10*95/ 

-1 

•175AQ.1 

/ •MFOmOaaG)  |R0  to  9 

f> 

CALL 

number 

1 0»90' 

-1 

.025AQ.1 

/ ElEM ( 2* J ) / n. « 3 1 

CALL 

SYMBOL 

1 1.00/ 

-1 

.175/0.1 

/ *MFD»aO,/3i  irO  to  q 

7 

CALL 

number 

( 1 • oc  • 

-1 

.025AQ.1 

/ ElEM ( 2# J ) / n. A 2 ) 

CALL 

Symbol 

1 1 . 05/ 

-1 

. 1 75a  0 . 1 

/ 'MFD'a0.a3i I go  to  9 

R 

CALL 

number 

1 1.00a 

-1 

•025*0,1 

/ELEH(2#jt/0./ 1 ) 

call 

symbol 

1 1.05/ 

-1 

• 1 75*0 . 1 

/ «MFO«  AO.A  3l 

9 CONTINUE 

return 

FnD 


I 

I 

r 


- AlO 


SuBS?9UTlNr  RC 


CBMMBN  X/Y/M.l 

.Q/T. J/ElEMj 1 6) 

X»X*0.75j 

CAlj 

PLB''’(X.Y,?i 

Y«Y-0.i5j 

CAl  1 

PLBT(X.Y,3l 

Y.Y*0.3i 

CALI 

PLflTfX.Y.?! 

x.x+0. 1 ; 

CA(  1 

PLBTf X,y.3l 

Y«Y»C.3i 

CALI 

PLBTf X.Y/2» 

Y.y-t.0. 1 

CALI 

PLBTf X.Y.31 

y»X*0.75j 

CALI 

PLBTfX.Y.Pi 

IF  ( L»LP«4  IELEM  (?*.i+7  ) f ?*  J + 7 ) *i  o • 0**A  • n 

Pifl  TB  ( li?#3«4i5/fti7#8iL 
1 CALL  NUMBERED. 75/0»37#0»liFLE'-uP*J-l  )<n.»-n 
CALL  SvMBftL(0»70/0»?0/0»l#  tPF'/0»>?)  I C,P  Tb  9 
? CALL  MJMBER»0.70>n.37#n.ljrLEM<?*J-l  wn..-n 

call  Symbol ( o»7c<n»BO#o • I *•  pf *# o» 1 rb  R 
3 Call  M)MBER(0.fe5i0»37#o*liFLEM(P#J-n#n..-1  > 
Call  SYMBBL(0.7C/n#!?0/0*I« 'PF'iOt#?)  I BB  Te  9 

A Call  MJ*^flFR(0.<i0in.37#o«l*FLEM(P*J-n/n.»-l  ) 
CALL  SYMBBL  f 0. 70/n.?0/0»  1 » iPF  ' /0»>?  > j BB  B 9 
5 call  MJM8ERr0.BC>0»37/n«1#FLEMjP*J-l 

Call  SV^RBLf0.ft'5iO»?0>0*l#  '^FD*/0»#3U  bp  tb  q 

b CALL  \UMBFR(0.BSin»37fO»1 iFLEM(P*j-t )in.»3» 

Call  symbbl ( o.ftRi n«20^ 0*  1 * t mfo * .0* #3 » 1 r,«  tb  q 

7 CALL  MJMBER(0.ftC^n.37/OtliFLEM(?*j.l )/n.»2) 
Call  BvMflPL  f O.B5/n»?C/0»  1 / •►^FO  ' iO»  / 3 ) ; r?P  TB  9 
« CALL  MJHBFRf  0*AC/0»37#o»1«F|.EM(?#j-i  wn.j  1 ) 
Call  svMBBL(0.f>R/n«20#o«i> 'mfo' »o«»3u  bb  tb  9 
9 CrNTImjE 

y.1.61  V-n.nj  CALL  PL6T(X»y#3) 

pfTurn 

FnD 


SuBRBUTiNr  9Cl 

CBM'^BN  X,Y«*^iL^TiT#J/El  EMMf  ) 
Imtegfp  Q 
Call  sc 


Y«Y»0.667?j 

CAl  1 

PLOT ( y, Y. ? 1 

Call  SLV(y.Y) 

Y.Y-0.667?! 

CAi.i 

PL8T(X,Y.2i 

y.y+2. 1 ; 

CALI 

PLBTf X.Y.31 

BB  TS 

1 Call  MJf^BFR  f 1 . 76/ -1  »02B#  0 • 1 > ELFM(  j ) *1  rjnO*  j 0«  « 1 I 
Call  SvMBBLf  1 tl7?5#n.1  / 'MH'/O./P)  iBB  tb  6 

? Call  number j 1 • 76/ - 1 *025/ 0 • 1 / ELEm ( 2# j ^ #1 nnO • # 0« # 1 ) 
Call  SyMBSL ( 1 .B6/-1 • 175/0. 1 / »MH» /O./?  ) iBB  TB  6 
3 Call  NliMBER(  1.76/-1. 025/0. 1 /ELEMf?,j)*innO./0..-n 
Call  symbol M .81 /-I . 1 75/C. 1 / 'MW i/O./P > iOB  tb  6 
A CALL  NljMBFRf  1.76/-1. 025/0. 1 /ELEMf2*jwn. .21 
CALL  Sy*^BBL  (1 .91  / -1 . 175/0.1  / 'H» /0«/ 1 1 iOB  Tb  6 

5 Call  number  ( 1 » 76/ - 1 . 025/ 0 . 1 . ELEM  ( 2#vJ  1 / 0 ..  1 1 
Call  symbol fi.9i/-i.i75/o.i/'w*/o./n 

6 CBNTINIIE 

return 

end 


- All 


SUB'^OUTINF  SLrRl 

JiEl.EMM  (S) 

I\TrGFR  OiT 
Call  si  cb? 


v.v-0 

CA|  1 

PLOT ( X, V, ? i 

x-X  + 0 

CAi  1 

PLOTf XjV.2l 

y.Y-0 

.I67?j 

CAi  1 

PLOT(X,v,?i 

Call 

Sl.V  ( Xi  V ) 

v.v-0 

.l67?j 

CAl  1 

PL0T(X.V,2i 

y .x-0 

CAl  1 

PLOTf X.ViPi 

V.  v*o 

.Bj 

CALI 

PL0T(X.V,2< 

x.x-0 

. 1 5 j 

CALI 

PLOTl X, V, 31 

x»x>0 

t3; 

CAl  1 

PL0T(X.Vi2l 

V.V4.Q 

• 1 ; 

CAl  1 

PL0T{X.V,3l 

x-x-o 

.31 

Cali 

PLOTf X.Yi2l 

X.XA’O 

.l5j 

CALI 

PL0T(X.Vi3l 

V.Y'AO 

• Bi 

Cali 

PLOT ( X«  Ya  2 1 

x.x^o 

• Pi 

CAl.l 

PL0T(X.Ya2i 

V.V-  1 

• i i 

Cali 

PL0T<XjYa3i 

v.v-0 

.B  J 

CAl  1 

PlOT<X#Ya2i 

V.V4? 

.1 ; 

CAli 

PL0TlX.Y.3i 

r,B  T0 

( 1a?a3/4,5J0 

1 

Call 

NU'^BER 

<?.36a-1 

. 02Ba  0.1 

*ELEM(2#j*gi*i OOO. 

A 0 . A 1 1 

CALL 

SV'^RBL 

f ? • 4 1 A - 1 

.175A0.1 

A 'mhi*o.a?i ior  to 

6 

? 

Call 

MJ^^BFR 

( ?.36a -1 

•025ao.1aElEM(2*J4Bi.iOOO» 

A 0.  A 1 1 

CALL 

SvRBBL 

f ?.46a-1 

. 175a  0. 1 

t 'MH«#0.a2) iRO  to 

6 

3 

call 

M/i^RER 

f ? • 36a  - 1 

. 02Ba  0.1 

«ElE^a(  2*j^0  1*1  OOO. 

A 0 . A - 1 

Call 

SvMRBL 

( ? • 4 1 A - 1 

•175A0.1 

. *MH»»0.a2) iRR  to 

6 

4» 

CALL 

MjMRER 

t 2 . 36  a - 1 

• 025  a 0.1  AELEM(2*j^Bwn»A2l 

CALL 

SvMRftL 

(P.BIa-1 

• 1 75 A 0 . 1 

A INI ,o»A 1 1 jOO  To  6 

B 

CALL 

NUMBER 

t 2 * 36a  - 1 

• 025 A 0.lAELEM(2*j*giin»Al  i 

CALL 

SvMrrl 

f 2.B1A-1 

. 175ao.  1 

A 'HI #0»A 1 1 

b 

IE  ( T. 

LF»4)ElEm(?*j) 

•ElEM(?* 

,) ) *1  0»0**6.O 

R0  Tfl 

( 7 J g # 

9a  1 riA  1 1 A 

12a 1 3a 1 4 1 T 

7 

CALL 

number 

( 1 . 45a  - 1 

.025ao. 1 

«ElEM(2.j).o,,-i  ) 

Call 

SVMR0L 

< 1 • 40a  - 1 

•175a0.1 

« 'PE*#0«a2> lOO  TO 

15 

A 

CALL 

MJMRER 

t 1 .40a  - 1 

•025a0. 1 

#ElEM(2*jiao.a*1  1 

CALL 

SvMBftL 

1 1 »40a -1 

•175A0.1 

A •PE«/C.a2) iRO  to 

1 5 

9 

CALL 

MJMRER 

f 1 .30a -1 

•025A0.1 

*ElEM(2*j)aO.»*1  1 

CALL 

SvMflflL 

( 1 . 35  A - 1 

.175a 0.1 

A *PE*a0.a21 iRO  TO 

15 

1 n 

CALL 

number 

< 1.2CA-1 

.025A0.1 

aElE^(2#j>ao..-1  1 

CALL 

symbol 

1 1.30A-1 

#175a0.1 

, »PE*/0.a?1  I ro  to 

15 

1 1 

Call 

number 

( 1 .OCa -1 

.025a0. 1 

.ElEM(2#ui#o.a4j 

Call 

SVMflflL 

1 1 . 1 5a  - 1 

• 175aq.  1 

$ *MED*A0.i3i ,rO  to 

15 

1? 

Call 

number 

f l«10A-l 

•025aq.1 

» ElE^ ( 2* J 1 A n. A 3 1 

CALL 

symbol 

f 1 .20a-1 

•175a0.1 

, *MEO' AO.A  3l irO  to 

15 

13 

call 

number 

( 1 .?0a-1 

•025ao. 1 

«ELEM(2«uiao.(2i 

CALL 

Symbol 

r 1 .?5a -1 

* 1 75a0. 1 

A *MEn»A0.A3i 1 GO  TO  15 

1 4 

Call 

number 

f 1 .PCa-1 

•025A0.1 

«ELEM(2#Ji/0.a 1 ) 

Call 

SYMBOL 

t 1 .?5a-1 

•175/0.1 

t *MEO*  AO.A  3l 

IR 

continue 

)CmS*Ot  Y»0*ni  CALL  PL»TfX#V#3) 

Return 

End 


) 


ff 
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SljBWRUTiNr  BLTR? 

X,  ,Q/T,  J/ELEMm  (S) 


ImTFGFR  G/T 
y»x+0.B'*7?j 

Call  si H(y,y ) 

X.X>*0»437?i 

V,V-0.i5j 

y.v*o.Bi 
y«x*o. 1 ; 
y»Y-0.3j 
y»v*o»i5i 
y »x-^0.B5  j 
GB  TP 

1 CALL  MjMRE;R(0.48»n 
Call  SyMBBL(0«83*0 
? Call  MJMBERfO»43iO 
CALL  Syf^BPLfOtBB/n 

3 Call  number ( o»48#n 
CALL  SyMBPLf 0#R3#0 

4 CALL  NuM0ER(O.48*O 
CALL  SyMB9L(0«93i0 

5 CALL  Nl)MBER(0.48#n 
CALL  SyMB0L(O.93/O 


If  ( L»LE»4  jElEM(  2*j+7  ) -ELEM  f ?*.j*7)#l  0.0*#A»n 


CA|  I PL0T(X,yi?^ 

CALI  PL9T<X,Y,2i 
CALI  PL9T(X,Y#3> 
cam  PL9T(X^Y,2i 
cam  PL9T(X,Yi3^ 

CALI  PL0T(X,Y>2< 

CALI  PL9T(X,Y,3» 

CAl  I PL9T(X,Yi2> 

14iO*1*EI  EM(?#J«1  )*1nnO««0»»I^ 
1A#0#1,  iG9  t9  f, 

lA<o.i,p-LEM(?*j-M*innO.*n.#l) 
14*0«1i iMH«iO»^?) iG9  70  A 
lAin«liELEM(2*J-n*inn0.iO«f»l) 
•MH'iO.i?) iGB  T0  h 
l4/ntl,FLEM(2*J.l  Mfl.iS) 
iM<  »Q.*  1 MG9  T9  6 
lAif)*!  #ELEM(2#J-1  mo.  » 1 ) 


r,9  T6 

(7#8*9* 10> 11 

/ 12/ 1 3/ 1 4 iL 

7 

CALL 

NUMBER ( 1 

• 0 

• 37/ 0 • 

1/ElEM!2*J47)/0^» 

•1  ) 

call 

SyMBBL  t 1 

• 50/  0 

• 18/0/ 

1/ tPE‘/0//?) 1 09 

Tp 

15 

8 

call 

NUMBER( 1 

• 50/0 

• 37/ 0 • 

1/FlEM!P#J*7|/0/» 

-1  ) 

call 

SyMB9L( 1 

• 50/0 

• 18/ 0 • 

1/ tpE»/0//Pl  1 09 

T9 

15 

9 

CALL 

NUMBER ( 1 

• 45/0 

• 37/  0 • 

1/ElEM!?»J47  W0^» 

-1  1 

Call 

SyMBBL ( 1 

• 50/0 

• 18/0/ 

1/ 'PE'/Oz/Pl 1 OP 

T9 

15 

1 n 

CALL 

Number  1 1 

• 47/0 

• 37  / 0 • 

1/El  EM!2«J^7)/0/» 

•1  > 

Call 

SYMB9L(  1 

• 50/0 

• 18/ 0 • 

1/»PE»/0//?M  09 

Tp 

15 

1 1 

CALL 

NUMBER(  1 

• 30/0 

• 37/0* 

l/ELEM!2*J47l/0/t 

^ ) 

call 

SYMB9L< 1 

• 45/0 

• 18/0^ 

1/  tMED»  /0»/3l  i 0»» 

TP 

15 

1? 

CALL 

NUMBER!  1 

• 35/0 

• 37/  0 • 

1 / el  Em ! ?*J^7 1 / 0^ . 

3| 

CALL 

SYMB9L  t 1 

• 45/0 

• 18/ 0 • 

1/ «MEO' /0»/3M  O'* 

T9 

15 

n 

Call 

NUMBER! 1 

• 40/  0 

• 37/0/ 

1/ElEM!?*J47W0/. 

2) 

call 

SyMBPL!  1 

• 45/0 

• 18/0/ 

1/ tMED* /0//3M  0« 

TP 

15 

14 

call 

NUMBER! 1 

• 40/0 

• 37/ 0 • 

l,ELEM!2*J+7i/0/. 

ll 

call 

SYMB9L! 1 

• 45/0 

• l8/0* 

1 / »mEO> /O*/ 3»  j 0« 

T9 

15 

IB  CBNTInuE 
X-2.0J  Y. 
RETURN 
END 


CALL  PL9T(X*Y#3> 


